Ultraslow quantum dynamics in a sub-Ohmic heat bath 
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We show that the low-frequency modes of a sub-Ohmic bosonic heat bath generate an effective 
dynamical asymmetry for an intrinsically symmetric quantum spin— 1/2. An initially fully polar- 
ized spin first decays towards a quasiequilibrium determined by the dynamical asymmetry, thereby 
showing coherent damped oscillations on the (fast) time scale of the spin splitting. On top of this, 
the dynamical asymmetry itself decays on an ultraslow time scale and vanishes asymptotically since 
the global equilibrium phase is symmetric. We quantitatively study the nature of the initial fast 
decay to the quasiequilibrium and discuss the features of ultraslow dynamics of the quasiequilibrium 
itself. The dynamical asymmetry is more pronounced for smaller values of the sub-Ohmic exponent 
and for lower temperatures, which emphasizes the quantum many-body nature of the effect. The 
symmetry breaking is related to the dynamic crossover between coherent and overdamped relax- 
ation of the spin polarization and is not connected to the localization quantum phase transition. 
In addition to this delocalized phase, we identify a novel phase which is characterized by damped 
coherent oscillations in the localized phase. This allows for a sketch of the zero-temperature phase 
diagram of the sub-Ohmic spin-boson model with four distinct phases. 

PACS numbers: 03.65.Yz, 64.70.Tg, 75.40.Gb, 05.30.Jp 



I. INTRODUCTION 

Open quantum dynamics^ describes relaxation and 
decoherence of quantum systems due to fluctuations in- 
duced by their environment. The noise characteristics 
determines the particular form of the relaxation dynam- 
ics and is specified by the spectral distribution J{oS) of 
the environmental modes. An important class is elec- 
tromagnetic noise which can generate simple frequency- 
independent damping. This is described by the well- 
known Ohmic environment J(u>) ~ w s with exponent 
s = 1. Even the simplest model of a quantum two- 
level system (TLS), the Ohmic spin-boson model, already 
shows nontrivial features, such as, e.g., a quantum phase 
transition of localization due to strong damping^. 

Recently, a different class of fluctuations has moved 
into the focus of attention^ - — , in which the low-frequency 
components are more strongly pronounced compared to 
the Ohmic case. This case is characterized by the sub- 
Ohmic distribution < s < 1, and the increased density 
of states of the low-frequency bath modes leads to pecu- 
liar and unexpected phenomena. The sub-Ohmic spin- 
boson model has been investigated by the flow equation 
technique^, which allows to calculate spectral proper- 
ties by equilibrium correlation functions. The results in- 
dicate a first order phase transition from derealization 
to localization at some finite value of the system-bath 
coupling^. Wilson's numerical renormalization group 
approach has allowed^ to identify a line of continuous 
boundary quantum phase transitions of the sub-Ohmic 
spin-boson model for all cases < s < 1. By the 
same method, weakly damped coherent oscillations on 
short time scales have been observed (but not further 



to be enhanced at the second order quantum phase tran- 
sition, which is indicated by a cusp in the entropy of 
entanglement 7 -. A continuous time cluster scheme based 
on Quantum Monte Carlo simulations has been used^ to 
determine the critical exponent at the quantum phase 
transition and to clarify a debate on the correctness of 
the quantum-to-classical mapping of the sub-Ohmic spin- 
boson model for s < 1/2. Recently, a novel numerical 
approach has been put forward^, which uses a sparse 
polynomial representation of the basis of the system-bath 
Hilbert space and which allows to obtain numerical re- 
sults based on exact diagonalization. These results were 
all obtained by numerical methods which do not rely on 
a pcrturbative approximation. On the other hand, these 
results also indicated^ that a perturbative approach-^— 
is by construction restricted to short times only, where 
it typically captures the renormalized coherent tunnel- 
ing, but has problems to describe the approach to ther- 
mal equilibrium in the long-time limit, see also Ref. [13 
for a related discussion. Finally, we mention the recent 
work of Wang ct alJ^ who discussed both the localization 
transition and the crossover between oscillatory and over- 
damped behavior. They find for all exponents < s < 1 
that with increasing coupling strength first the oscilla- 
tory behavior becomes overdamped and only then (for 
larger couplings) the transition to localization follows. 

Sub-Ohmic fluctuations play a central role for the 
quantum behavior of nanomechanical devices^ due to 
the coupling of flexural modes to many surface defects 
of the resonator. Moreover, an effective sub-Ohmic spin- 
boson model is studied as the archetype local quantum 
critical system, showing a quantum critical transition 
(QCT)±&. The nature of the QCT in heavy Fermion 
metals is crucial to understand quantum phase transi- 
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An important experimental puzzle 



turn entanglement of the spin and the bath has shown is whether a non-linear scaling of the dynamic spin sus- 
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ccptibility x(w) ~ (u>/T) v with v < 1 occurs 1 ^ (here, w is 
the frequency), which allows for the experimental iden- 
tification of the QCT in the spin density wave. Follow- 
ing a five-dimensional 4 field theory, the corresponding 
quantum critical point in three dimensions can be iden- 
tified by a Gaussian fixed point. A violation of the linear 
scaling v = 1 would indicate an interacting fixed point 
with spin damping being present, while a nonlinear scal- 
ing v < 1 follows from assuming a non-interacting fixed 
point (with spin damping being absent). The question 
of the deviation from linear scaling is at present not fi- 
nally answered 1 ^. Moreover, the sub-Ohmic spin-boson 
model applies to superconducting nanoscale device a 20 ' 21 , 
such as qubits in contact with RLC transmission \mesM. 
Finally, the limit of s — > corresponds to the important 
class of f //— nois o 20 ' 23 i 24 . Often 1//— noise is generated 
by two level fluctuators as known from glasses^ where 
they cause ultra slow (aging) dynamics 2 ^—. Due to the 
dominant low-frequency modes, the dynamics is notori- 
ously difficult to determine, since strong non-Markovian 
effects even for relatively weak coupling arise. 

In this work, we show that a sub-Ohmic heat bath 
generates an effective dynamical asymmetry for an in- 
trinsically symmetric quantum TLS (spin— 1/2), which 
leads to a transient shifted quasiequilibrium for the TLS, 
which then itself decays extremely slowly. This tran- 
sient symmetry breaking occurs in the delocalizcd phase 
of the TLS, which follows the slowly changing quasiequi- 
librium. Moreover, it depends on the precise form of 
the initial condition and only occurs for the case when 
the spin— 1/2 is initially shifted and the bath is initially 
thermalized to the shifted spin. We establish a simple 
qualitative picture based on a perturbative analysis in 
terms of an asymmetric spin-boson system in order to 
interpret the dynamics of the spin polarization at short 
times. There, the dynamical asymmetry shows an initial 
fast exponential decay towards a quasiequilibrium with a 
finite quasistatic asymmetry. Since the system is globally 
symmetric, a subsequent decay towards global thermal 
symmetric equilibrium follows in which the dynamical 
asymmetry vanishes as expected, as long as the system 
is in the delocalizcd phase. Our numerical data allow us 
to conclude that the time scale of this final decay is much 
larger than any simple time scale, such as given by the 
inverse spin splitting A -1 , the inverse temperature T _1 
or the inverse cut-off frequency w" 1 of the bath. More- 
over, we find that it is much larger than the observed 
decoherence time 7^~ 1 . The asymptotic regime cannot 
be described anymore in terms of a simple perturbative 
picture. 

The feature of the ultraslow relaxation dynamics to- 
gether with the fast initial decay towards a quasiequi- 
librium resembles qualitatively the dynamical response 
of coarsening or a spinglass-like dynamics^ to some ex- 
tent. We should, however, emphasize that the numerical 
scheme is limited to time ranges which are too short to 
study the full asymptotic nature of the ultraslow decay, 
apart from stating its existence and some central fea- 



tures, such as its dependence on the model parameters. 
A full study of quantum critical features lies beyond the 
scope the approach. In general, the ultraslowly decaying 
transient asymmetry is shown to be more pronounced for 
stronger system-bath coupling and smaller exponents s, 
which clearly illustrates the quantum many-body nature 
of the dynamical asymmetry. In addition to this, we find 
coherent oscillations in the noncquilibrium polarization 
of the TLS in the localized phase where in equilibrium 
all coherence is suppressed. These findings allow us to 
sketch the zero-temperature phase diagram with four dis- 
tinct phases. We argue that our results may have imme- 
diate experimental implications, since the slow transient 
asymmetry will dominate all spin fluctuations and thus 
the true equilibrium spin fluctuations, which are relevant 
for quantum critical behavior and which require asymp- 
totically long times, may be difficult to access experi- 
mentally. In turn, the slow dynamical asymmetry should 
also be taken into account by quantum engineering of 
nano devices. 

To structure the paper, we introduce the model in 
the next section and then discuss briefly the necessary 
technical modifications of the numerical exact quasi adi- 
abatic propagator path-integral to tackle the spin dynam- 
ics with factorized and shifted or polarized bath initial 
conditions. In Section IIV1 we develop our simple quali- 
tative picture which quantitatively describes the forma- 
tion of the dynamically generated asymmetry for the spin 
at intermediate times. In Sec. [Vj we sketch a complete 
phase diagram of the sub-Ohmic spin-boson model which 
summarizes our findings and which includes four distinct 
phases, before we finish with the conclusions. 



II. THE MODEL 

The spin-boson Hamiltonian is given by (ft = 1), 
H = Hs + Hsb + Hb with H$ = \ Acr^ — \^o~z with tun- 
neling element A and asymmetry e. Moreover, Hsb = 
~lf Eft ^kQk and H B = \ Efc (p* + ^M) , wh erc q k and 
Pk are the position and momentum operators of the bath 
normal mode with wave vector k. The bath is completely 
described by the spectral function 

J(w) - Y^ 6 ^ - - 2au>l~'w'e~ u/Ue (1) 
ft k 

with the bath cut-off frequency uj c and the dimcnsionlcss 
coupling a. In this work, we set cu c = lOA^i. The time 
dependent polarization P(t) := (a z )(t), after initial full 
polarization, P(0) = 1, shows coherent damped oscilla- 
tions and for an asymmetric two-level system addition- 
ally an exponential decay. In the Ohmic case, s = 1, the 
model displays a zero temperature quantum phase tran- 
sition at a c = 1 between a delocalizcd phase, where the 
spin can tunnel between its two states, and a localized 
phase. Even in the delocalized phase at a < a c , co- 
herent dynamics is found only at weak coupling a < 5. 
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FIG. 1: Polarization P(t) for T = 0.1A, s = 0.25 and u c 
10A for factorized and polarized bath initial conditions. 



For | < a < 1, ovcrdamped relaxation occurs with a 
dynamic crossover at a = ad — \- In the sub-Ohmic 
regime, both a continuous (in s) localization transition 
and dynamic crossover are known 3 ^^^, but the inter- 
play of both over the full range of exponents s has not 
been discussed and its phase diagram is unidentified. 

To determine P(t) , we use the numerically ex- 
act real-time quasiadiabatic propagator path integral 
(QUAPI) 32 ' 33 . We are numerically restricted to finite 
times and cannot investigate the localization quantum 
phase transition, but rather the dynamic crossover be- 
tween oscillatory and overdamped dynamics. We calcu- 
late P(t) for a symmetric TLS (e = 0) in the delocalized 
phase for three different exponents s = 0.25, 0.5 and 0.75 
and focus mainly on low temperatures. We study the cou- 
pling range from weak to strong (delocalized/localized 
phase) as long as the dynamics is oscillatory. 



III. QUAPI FOR THE SHIFTED BATH INITIAL 
CONDITION 



The dynamics of the sub-Ohmic spin-boson model is 
rather sensitive on the way how the system and bath 
are initially prepared, in particular with respect to their 
mutual interaction. It enters the formalism when the 
bath degrees of freedom are averaged over. To under- 
stand this, consider an experiment where a symmetric 
spin is in thermal equilibrium with the bath at temper- 
ature T and at initial time t = 0, the spin density op- 
erator ps is prepared such that (o~ z )(t = 0) = I. Then, 
the bath is in thermal equilibrium to (o~ z )(t' < 0) = 0, 
resulting in the factorizing initial condition for the to- 
tal density operator at initial time, po = Ps ® Pbo 
with pbo = exp(— Hb/T)/Zb with ks = 1,-Hb = 
iEfcCPfc + WfcOfc) and Z b = Tr cxp(—Hb/T). A dif- 
ferent case is to start preparing the spin long time be- 
fore the experiment starts in order for the environment 
to thermalize to the shifted spin. This results in the 



shifted (correlated) initial condition pg = ps <8> psq with 
PBq = cxp(-{ff B + H SB \ az=1 )/T}/Z Bq , where H SB = 
-it T. k ^q k and Z Bq = Trexp(-{H B + Hs B \„ M =i)/T}. 
All recent studies 3 -^^ were based on the shifted initial 
condition. 

The coupling between the spin and a single bath mode 
scales with 1 / y/N, with TV being the number of modes 
and thus is infinitely small for a continuum of modes. 
Therefore the difference between pso and ps q is small 
and believed to have little influence on the spin dynam- 
ics itself. This assumption is at the heart of open quan- 
tum dynamics since only an environment which is large 
enough to stay unaffected by a small quantum system can 
be averaged over without loss of crucial physics. Lucke 
et alM: showed that for an Ohmic bath, both initial con- 
ditions lead to sizeable differences for the spin dynamics 
only on very short times of the order of l/u c . Thus, 
the spin dynamics which lives on a time scale of 1/E, is 
unaltered as long as E <C uj c . 

The relevant parameter to distinguish both cases is 
the energy difference between the two initial conditions 
of the bath. It is estimated by the reorganization energy 
Er = J Q °° dui^^- and is dominated by the low frequency 
modes. As they are weakly present in the Ohmic case, 
their influence is negligible when A -C w c . In contrast, for 
a sub-Ohmic bath, any reorganization is slow, implying 
that the bath initial condition has a crucial influence on 
all time scales of the system dynamics. 

In order to treat this situation properly, we adopt the 
QUAPI scheme, which has originally been formulated for 
factorizing initial conditions^, also to correlated initial 
conditions. QUAPI is an efficient, iterative and numeri- 
cally exact method, in which the real-time quantum evo- 
lution operator is sliced into Af Trotter slices of duration 
At. For the factorizing initial condition, the bath influ- 
ence is captured by influence functional 

{M k \ 
fe=0fc'=0 J 

with the bath correlators rjkk' defined in Ref. |32| and the 
constant spin coordinate Sk = \o z = ±| at the time in- 
terval [(k — ^)At, (k + i)At]. At any finite temperature, 
all bath correlations decay exponentially fast at asymp- 
totically long times 3 ^, which allows for a truncation of the 
bath-induced memory time beyond a certain time span. 
This corresponds in the above sum to setting rjkk' = 
when \k — k'\ > k max for a given memory length k max . 
Finally, convergence with respect to increasing k rnax has 
to be ensured. Inside the finite size memory time win- 
dow, all correlations are taken into account exactly and 
the corresponding path sum is carried out deterministi- 
cally. 

For the shifted bath initial condition, we follow Ref. 
[35l to modify the influence functional I q = iocxp(— $ 9 ) 
with 
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FIG. 2: Dynamically generated asymmetry e(t) for T = 0.1A, 
s = 0.5 and ui c = 10A for the polarized bath initial condition. 
The exponential decay is clearly visible by the semilogarith- 
mic scale. Inset: Parameters extracted from a fit to Eq. ([2]) 
as a function of a (symbols); solid lines are linear fits. 



with F(T,w,t') = J 1/T rfrsinh[w/(2T) + i(t' + ir)] = 
2 sinh[w/(2T)] cos(ut')/ui, the spin s(i) = ^cr z (t) and the 
shift q 

■ sin(wi) 



i for the shifted initial condition. After Trotter 
slicing, we obtain with = Jn^i^ t duj — 



The bath correlators ry^/ are all two-time correlation 
functions between the times k 1 At and kAt. In contrast, 
Vq,k describes a correlation between time kAt and the 
initial time t = 0. Within the QUAPI scheme, neglect- 
ing all correlators with \k — k'\ > k max would imply to 
also neglect rj qt k at times t > k max At and no sizeable 
difference for the two initial conditions would be found. 
However, since rj q ^ depends only on a single running time 
index k, we can incorporate it completely by calculating 
the necessary r\ q ^ in each time step during the iteration. 
Thus, the shifted bath correlations are taken fully into 
account, independent of k ma x- 

The two initial preparations induce qualitatively dif- 
ferent behaviors of the dynamics of the sub-Ohmic spin- 
boson model. Fig. [T] shows the polarization Pit) for 
the factorizing initial condition for a symmetric spin 
(e = 0) for different coupling strengths for T = 0.1A 
and s = 0.25 (a c = 0.022^). We find damped oscillations 
where damping grows with increasing a. As expected 
from a Bloch-like picture, the oscillations occur around 
zero average. Surprisingly even in the localized phase, 
P(t) shows quantum coherent oscillations out of equi- 
librium. In the localized phase the bath modes are too 
slow to localize the out-of-equilibrium TLS fast enough 
and thus damped coherent oscillations still occur at short 
times. 

In contrast to that, Fig[T]also shows P{t) for the shifted 
bath initial condition. As before, we find damped oscilla- 
tions with a damping constant increasing with a. How- 
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FIG. 3: Initial (e(t = 0) = e ) and final (e(t = 50/A)) dy- 
namical asymmetry, tunneling matrix element Aat and deco- 
herence rate 72, At at T = 0.1A for the polarized bath initial 
condition for s = 0.25 (a), 0.5 (b) and 0.75 (c). Note that 
in (c) the data (magenta squares) for e(t = 50/A) have been 
scaled up to increase visibility. 



ever, in striking qualitative difference, the oscillations 
occur around a finite value P av > 0, even for coupling 
strengths well in the delocalized phase, which indicates 
an intermittent symmetry breaking. Note that neither 
P av > nor the long-time value P(t = 25/ A) nor the oc- 
currence of oscillations allow to distinguish between the 
delocalized and the localized phase. The latter results are 
in line with those of Ref . @. The dependence on the initial 
preparation, however, has not been disclosed before. 



IV. ULTRASLOW DECAY OF THE EFFECTIVE 
ASYMMETRY 

In order to study the effect of the shifted oscillations 
more quantitatively, they can be formalized in terms of 
an effective time-dependent asymmetry e(t) for the TLS 
which generates a finite polarization P av . The polariza- 
tion of a TLS with a static asymmetry e at low temper- 
atures and weak coupling is given^ by 



P(t) 



{cos(Et) + -^sm(Et)} 



,-71* 



72 

E 
e 
E 



-72* 



(l-e^ 1 *) 



(2) 



with E = \J A 2 + ( 2 . Certainly, for our intrinsically sym- 
metric case, at asymptotically long times, the asymmetry 
has to vanish, while our data reveal that the asymmetry 
is generated at short to intermediate times, rendering e 
effectively time-dependent. Nevertheless, we expect that 
particularly at short times, the dynamics is well described 
by a lowest order perturbative description in the tunnel- 
ing strength, yielding an initial exponential decay also 
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for the bath correlations. Thus, it is suggestive to use 
the heuristic approach e — > e(t) = e + e i cx p(~ 7b f° r 
the dynamical asymmetry and E = E(e(t)) in Eq. (j2J|. 
Clearly, this approach breaks down at asymptotically 
long times since it would yield a finite static asymptotic 
bias e{t — > oo) = eo, which cannot occur in a globally 
symmetric system. Nevertheless, we use this heuristic 
approach to show the features of the dynamical asymme- 
try on a qualitative level. 

Fitting our data for P(t) at short to intermediate times 
with Eq. ([2]) allows us to extract eo,ei and 7b- Fig. [2] 
shows the result of the fit for e(t) for T = 0.1A, for 
different a and s = 0.5 with critical coupling for local- 
ization a c ~ 0.1065^. The results of the fits with the 
effective dynamical asymmetry are perfect (not shown) 
on the considered time span up to t = 50. This approach 
will no longer hold for truly asymptotic times t — > 00. 
With increasing a, the dynamical asymmetry (both at 
t = and t = 50A -1 ) increases, but its short time decay 
rate jb decreases slightly. The inset of Fig. [2] shows the 
extracted parameters vs. a and linear fits to the data. 

A simple qualitative picture can be obtained by assum- 
ing that the initially shifted bath modes generate this 
transient dynamical asymmetry. The slow modes drag 
the TLS to its original position, similar to the cage ef- 
fect in the fractional Langevin equation 3 ^. This causes 
a non-zero effective polarization P av 7^ and a finite 
asymmetry both of which are clearly dynamically gener- 
ated by the sub-Ohmic bath correlations. In general, the 
symmetry breaking is enhanced with increasing coupling, 
emphasizing the many-body character of the asymmetry. 
At asymptotically long times, we expect the bath to re- 
lax slowly towards global equilibrium with the TLS be- 
ing in either of its two symmetric eigenstates in which 
(&z)eq = and the asymmetry will thus slowly vanish. 
The decay of the initial coherent oscillations shown in 
Fig. [1] occurs on a time scale of the same order as the 
period of the oscillation, but is very fast compared to the 
decay of P av which we actually cannot fully observe but 
clearly conjecture 3 - 7 -. This effect depends on the initial 
preparation of the bath and, moreover, it is absent in 
an Ohmic bath^i. We note that the feature of an ini- 
tial decay with a time scale of the order of a microscopic 
decay rate of the spin and a subsequent ultraslow de- 
cay towards a global equilibrium resembles that of coars- 
ening or of spin glasses 3 ^ after a quench 3 ^. Although 
the sub-Ohmic model has no explicit spatial disorder, its 
low-frequency noise implicitly mimics disorder, similar to 
l//-noise. The latter results from localized fluctuators 
with broadly randomly distributed parameters^ 3 -. The 
slow transients might render the sub-Ohmic spin-boson 
model a toy model to study quantum glasses^. 

Fig. [3] presents results for the fit of the initial dynam- 
ical asymmetry e(t = 0) = eo + ei, the asymmetry e(t = 
50/ A) for the maximal simulation time, and the tunnel- 
ing clement Agt and dccohcrence rate 72. fit vs. a/a c for 
s = 0.25,0.5 and s = 0.75, with a c ~ 0.022,0.1065 and 
a c ~ 0.3^, respectively. Note that we have set 71 = 272, 




temperature T [A] 



FIG. 4: Main: Ratio of tunneling element A/« and decoher- 
ence rate 72, fit vs. T for a/a c — 0.5. Inset: The correspond- 
ing final asymmetry e(t = 50/ A) = e j it +eij it exp(--y B ,fitt) 
vs. temperature. 



which is strictly valid for a symmetric system^ and that 
we have ignored the additional correction term for the 
rate due to a finite asymmetry. All data show the initial 
oscillatory behavior except for s = 0.75 and a = 0.7a c 
which exhibits overdamped dynamics. For s = 0.25, the 
asymmetry is smoothly increasing for larger a up to the 
phase transition and even beyond. At the same time, 
the tunneling matrix element is hardly influenced and 
larger than the decoherence rate which increases with 
coupling. Thus, the presence (absence) of oscillations is 
not a sufficient condition for the existence of the delo- 
calized (localized) phase. For s = 0.5, the asymmetry is 
generally smaller, but also increases with increasing a. 
For a > 0.9a c , the decoherence rate exceeds the tun- 
neling element hinting at the crossover to overdamped 
behavior for stronger couplings. For s = 0.75, the dy- 
namical asymmetry is small, depends nonmonotonously 
on a and vanishes for a > 0.7a c . The tunneling ma- 
trix element also vanishes there, indicating the crossover 
to overdamped dynamics. Thus both the overdamped 
dynamics and the dynamical asymmetry are closely con- 
nected. 

Thus, we can summarize our physical picture in the 
following way: In general, when the TLS is fully polar- 
ized and released at t = 0, it starts to oscillate coher- 
ently and resonant bath modes damp the oscillations. 
Since the slow bath modes are not yet in equilibrium 
with the evolving TLS, but still shifted, they cause a 
drag for the TLS to its initial state, resulting in the 
formation of a dynamical asymmetry. This decays to 
zero on an ultraslow time scale since the global equi- 
librium is symmetric. Thus, the local quantities decay 
first to a local quasicquilibrium generated by the imme- 
diate environment. This oscillatory decay at short to 
intermediate times is well characterized by exponential 
decay rates 71/2 (see Eq. ([2]), data not shown explicitly) 
and a dynamical asymmetry e(t), which itself varies with 
time. Following our simple perturbative line of reason- 
ing, e(t) = eo + £i exp(— 7b t) consists of two parts. One 
part, ei exp(— 7si), decays exponentially with rate jb 
and quickly vanishes as time increases. The second part 
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spectral exponent s 

FIG. 5: Sketch of a phase diagram of the sub-Ohmic spin- 
boson model at T — 0, see text. 

eo decreases only very weakly with growing time which 
is revealed by our numerical data, but has to vanish in 
the asymptotic limit t — > oo due to global symmetry. 
Hence, only at very large times, a global equilibrium 
state emerges. At stronger coupling, strong correlations 
between the TLS and resonant modes generate an over- 
damped dynamics. The responsible modes are not the 
slow modes, thus the dragging effect is suppressed. 



V. DYNAMIC CROSS-OVER AND SKETCH OF 
A PHASE DIAGRAM 

With increasing temperature, the dynamics becomes 
overdamped as well. Fig. @] shows the ratio Afi t /72,fit for 
a/a c = 0.5. In general, it decreases with increasing T. 
Only for s = 0.25 at T = 3A, we observe overdamped 
behavior. The data shown in Fig. @] indicate that the 
crossover temperature at which the overdamped regime 
start increases with decreasing exponent s (at fixed rela- 
tive coupling a/a c = 0.5). The final asymmetry (inset of 
Fig. 2]) decreases with temperature in a monotonic way, 
also around the crossover temperature. 

The findings can be summarized in a phase diagram 
sketched in FigJS] (at T = 0). The black circles mark the 
equilibrium localization phase transition known in the 
literature^ (note that the numbers refer to the values of 
a along the ordinate, the scale, however, is not linear). 
The second solid line with the squares mark the dynamic 
crossover from coherent to overdamped dynamic s 11 ' 40 . 
The delocalized (weak-coupling or symmetric) phase is 
reflected in the coherent oscillations of P(i), in which we 
have unexpectedly found a transient symmetry breaking 
for all s < 1. For l/2<s<l,a dynamical crossover to 
the overdamped delocalized phase occurs without tran- 
sient symmetry breaking, before the localized phase is 
reached for large a. Unexpectedly, for s < 1/2, no over- 



damped dynamics occurs but the nonequilibrium polar- 
ization still shows quantum coherent (damped) oscilla- 
tions in the localized phase, where in equilibrium coher- 
ence is completely suppressed. Although we could not 
obtain converged results for stronger couplings, we ex- 
pect that the oscillatory dynamics deeper in the localized 
phase should vanish. 



The phase diagram summarizes our two major results: 
(i) In the delocalized phase where the dynamics is oscil- 
latory, the initial preparation of the bath strongly influ- 
ences the dynamics. Thus, the dynamics in this phase 
is dominated by the ultra-slow sub-Ohmic bath modes. 
Any experiment addressing this phase will have to deal 
with these unusual long-time transients and true equi- 
librium quantities will be hard to access. Accordingly, 
one cannot expect to observe the spectrum of the equi- 
librium fluctuations when approaching the localization 
phase transition from this phase. This has to be taken 
into account when analyzing the data for the dynamic 
structure factor of heavy fermion systems within a pic- 
ture of local quantum criticality. The latter assumes the 
local critical fluctuations to be of sub-Ohmic nature, (ii) 
Moreover, the phase denoted as "localized-oscillatory" 
in Fig. [5] is characterized by an oscillating polarization, 
although the spin is in the localized phase where in equi- 
librium any oscillation is suppressed. This phase is novel 
and, to the best of our knowledge, has not been revealed 
or discussed before. Note that it is not an equilibrium 
phase. Strictly speaking, it is just a part of the equilib- 
rium phase denoted as "localized" , but from the view- 
point of the real-time dynamics, the two regions (one 
with oscillatory dynamics and one without) are distinct. 
It is known that for s > 1/2 no oscillatory dynamics oc- 
curs for couplings stronger than the critical coupling of 
the crossover to overdamped dynamics. Accordingly, in 
the localized phase (which occurs here at even stronger 
couplings), no oscillatory dynamics is observed. We be- 
lieve that even for s < 1/2, the oscillatory behavior will 
be suppressed at very large couplings. Hence, we conjec- 
ture the dashed phase separation line between localized- 
oscillatory and the localized-overdampcd region in Fig. 

El 



We note that experimental studies in this phase might 
observe a dynamical behavior which does not allow to dis- 
tinguish it from behavior found in the delocalized phase 
unless the experimental time scale exceed the time scale 
of the ultraslow decay in the delocalized phase since in 
the localized phase the asymmetry will not decay at any 
time. Although the numerical scheme does not allow to 
observe the ultraslow decay, we can estimate it to be at 
least two orders of magnitude slower than the intrinsic 
time scale A -1 of the spin as we observed finite asym- 
metries up to times 50 A" 1 (see Fig. [2]). 
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VI. SUMMARY 

To summarize, we have found a distinct ultraslow tran- 
sient dynamics of the polarization in the sub-Ohmic spin- 
boson model. A dynamically correlation-induced asym- 
metry is generated, which mimics a local quasicquilib- 
rium and generates a symmetry broken phase at inter- 
mediate times. It occurs due to the dominating low- 
frequency bath modes and only for an initially polar- 
ized thermal bath. The latter refers to the case when 
the system is initially prepared in a non-thermal polar- 
ized state and the bath can thermalize to that before 
relaxation to the global thermal equilibrium can start. 
Technically, this kind of initial preparation can be in- 
cluded by adopting the numerical quasiadiabatic propa- 
gator path integral scheme. For this, an additional term 
in the Fey nman- Vernon influence phase is included which 
only depends on the initial and the momentary observa- 
tion time, but not on the times in between. Hence, the 
general strategy of a finite memory time in QUAPI is not 
affected and this term can be taken into account in an 
exact manner. With that, we can show that the dynamic 
crossover between oscillatory and overdamped behavior 
is connected to the transient symmetry breaking and not 



to the localization quantum phase transition. 

Finally, we comment on an important possible conse- 
quence of our findings for identifying the quantum crit- 
ical transition in heavy Fermion system o 17 ' 18 . An effec- 
tive sub-Ohmic spin-boson model is used for thisi£. At 
present , the non-linear scaling of the dynamic suscep- 
tibility ~ ( UJ /T) l/ with v < 1 is used as a crite- 
rion^. However, our results suggest that the ultraslow 
relaxation dynamics in the sub-Ohmic spin-boson model 
requires to carefully measure the relaxation dynamics at 
very long times. In fact, this could render an experimen- 
tal observation of the equilibrium localization quantum 
phase transition difficult in general. For the same rea- 
son, the nonlinear scaling of the dynamical susceptibility 
has to be evaluated very carefully in order to ensure the 
existence of local quantum criticality in heavy Fermion 
systems. 
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